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Abstract 

Secluded dark matter models, in which WIMPs annihilate first into metastable mediators, 
can present novel indirect detection signatures in the form of gamma rays and fluxes of 
charged particles arriving from directions correlated with the centers of large astrophysical 
bodies within the solar system, such as the Sun and larger planets. This naturally occurs if 
the mean free path of the mediator is in excess of the solar (or planetary) radius. We show 
that existing constraints from water Cerenkov detectors already provide a novel probe of the 
parameter space of these models, complementary to other sources, with significant scope for 
future improvement from high angular resolution gamma-ray telescopes such as Fermi-LAT. 
Fluxes of charged particles produced in mediator decays are also capable of contributing a 
significant solar system component to the spectrum of energetic electrons and positrons, a 
possibility which can be tested with the directional and timing information of PAMELA and 
Fermi. 
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1. Introduction 



The search for weakly interacting massive particles (WIMPs) as a component of non-baryonic 
dark matter has become a focal point of modern particle physics [1]. There are several 
complementary experimental and observational approaches to WIMP detection [2]. Direct 
detection experiments probe the terrestrial scattering of WIMPs with nuclei and typically 
require low radiation environments to keep backgrounds under control. High energy colliders 
such as the Tevatron and the LHC offer the possibility of producing WIMPs and measuring 
their properties in the laboratory, provided the challenging missing energy signatures can 
be disentangled. Indirect searches for dark matter annihilating into gamma and cosmic rays 
in the galactic halo are also promising, although susceptible to various, often uncertain, 
astrophysical backgrounds. Finally, neutrino telescopes such as Super-Kamiokande and Ice 
Cube can search for indirect evidence of the annihilation of WIMPs captured in the core of 
the Sun and the Earth, in the form of an observable muon signature arising from neutrino 
charged current scattering in the detector. While the latter two examples are well-known 
indirect signatures for any thermal relic WIMP dark matter candidate, more generic WIMPs 
forming part of a larger dark sector can lead to further novel signatures. In this paper, we 
demonstrate that models of secluded dark matter |3] present an additional observational 
possibility: high-energy gamma rays and charged particles arriving from a direction tightly 
correlated with the centers of the Sun, Earth and other planets. Such novel signatures can 
be effectively probed with the powerful new generation of gamma ray telescopes. 

The primary feature of secluded models of dark matter [3J is a two-stage dark matter 
annihilation process: WIMPs annihilate first into metastable mediators, which subsequently 
decay into Standard Model (SM) states. This breaks the more-or-less rigid link between 
the size of the WIMP annihilation and WIMP-nucleus scattering cross sections. It has been 
shown that a small mass for the mediator allows for new phenomenological possibilities in 
the form of enhanced WIMP annihilation at small velocities [U [S] that may help to explain 
various astrophysical anomalies, e.g. the positron excess observed by PAMELA above 10 
GeV [6j and perhaps the unexpectedly hard electron spectrum observed by Fermi above a 
few hundred GeV |7]. Furthermore, a relatively small mediator mass kinematically removes 
heavy SM particles from the final state ^ [5], reconciling these effects with the absence 
of any enhancement in the cosmic ray anti-proton signal [8]. The lifetime of the mediator 
is essentially a free parameter, limited only by the Big Bang Nucleosynthesis bounds of 
r < 1 s. If this lifetime is rather long, the decay of the mediator will occur a long distance 
away from the point of the original WIMP annihilation. Denoting the WIMP particle x ^i-nd 
the mediator particle V , assuming XX ~^ 2V^ as the main annihilation channel, and taking 
my <^ m^, we arrive at the following estimate for the mediator travel distance: 

L = crv^7v^ = 3xl0^kmx^x^. (1) 

Such large boosts 7^/ = m-^/my are easily achieved if the the dark matter mass is near the 
electroweak scale and the mediator mass is below a GeV. With regard to the annihilation 
of WIMPs captured within the Sun, one can see that this distance may very well exceed 
the solar radius {Rq = 6.96 x 10^ km) in which case, unlike conventional WIMPs, most of 
the decay products will not be absorbed. For somewhat shorter lifetimes, the interesting 
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Figure 1 : A schematic illustration of the new indirect detection signature of secluded WIMPs captured in 
the solar core, annihilating to metastable mediators and leading to an electromagnetic flux: 7, e^, /i^, • • • . 
Sensitivity to conventional WIMPs arises only through annihilation to neutrinos. 

possibility emerges of mediators produced by WIMP annihilation in the center of the Earth 
decaying directly into charged particles within neutrino telescopes. A schematic illustration 
of this new indirect mechanism of probing secluded WIMP models is shown in Fig. [H 

A sub-GeV mass mediator will decay predominantly into light states, such as pions, 
muons, electrons, gammas and neutrinos. As they are produced in the decays of highly 
boosted mediators, these decay products will be tightly correlated with the direction to the 
original WIMP annihilation point, with a typical angular size 9y ~ l/7y. While for charged 
particles this correlation can be reduced by the magnetic fields encountered on the way to 
the detector, the directionality of gammas and neutrinos is unaltered. Modern gamma ray 
telescopes enjoy an angular resolution much better than a degree, which may be exploited to 
enhance the gamma ray signal-to-background. In particular, a notable source of background 
to this signature is the generation of gamma rays via cosmic rays impinging on the Sun (see, 
e.g. [HI [in])- Such background gammas will typically display a much softer spectrum than 
gamma rays from mediator decays and will have no specific correlation with the solar center 
where most of the secluded WIMP annihilation is expected to take place. The main gain in 
sensitivity in detecting gammas from the decays of metastable mediators, as compared to a 
more conventional search for a highly-energetic neutrino signal, may come from the increase 
in efficiency. While the detection of multi-GeV neutrinos requires their conversion to muons, 
which means a loss in efficiency of around ten orders of magnitude, the efficiency of detecting 
gammas created outside of the solar radius can be order one. 

In addition to gamma rays, many secluded WIMP models with relatively light mediators 
are destined to produce a significant fraction of leptons in the final state, and thus we are 
naturally led to the question of whether mediator decays outside the solar /planetary radii 
are capable of contributing significantly to the fiuxes of electrons and positrons seen by 
PAMELA and Fermi. It is tempting to pursue the notion that these anomalies may in fact 
have a local origin in the solar system, powered by the annihilation of secluded dark matter 
trapped within solar bodies. This is an intriguing possibility, as it would offer a new dark 
matter interpretation of these signatures that does not rely on galactic WIMP annihilation, 
which generally requires a significant boost factor in the annihilation cross section. In this 
context, a less extreme hierarchy between the WIMP and mediator mass may be feasible. 
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and this idea is akin to using local substructure to enhance the charged particle flux, albeit 
with a source (the Sun) which is extremely local and well-understood. These electrons and 
positrons, like the gamma rays discussed above, should be strongly correlated with the center 
of the Sun, particularly in the high energy range where the effect of the magnetic fields of 
the Sun and Earth is less significant. Such a hypothesis appears straightforwardly falsifiable 
using directional and timing data from PAMELA and Fermi. Finally, if charged particles 
are the primary signatiuc, this immediately leads to a minimal gamma ray flux from the 
associated 3-body decays that is also directionally correlated with the center of the Sun, and 
thus testable by Fermi-LAT. 

In this paper we analyze the feasibility of detecting electromagnetic particles, 7, e^, • • • , 
arising from the delayed decays of metastable mediators produced through WIMP annihi- 
lation in the deep interior of the Sun and other planetary bodies within the solar system. 
Our primary focus is on the capability of modern gamma ray telescopes to search for these 
local annihilation signatures and the possibility of distinguishing the signal gammas from 
the solar backgrounds produced by cosmic rays. We demonstrate that in certain corners of 
the secluded dark matter parameter space, the gamma-ray signature from the center of the 
Sun is potentially the most sensitive probe, superior to direct detection and other indirect 
signatures. This is especially the case for models in which WIMP-nucleus scattering is dom- 
inated either by spin- dependent interactions or proceeds through an inelastic transition to a 
nearby excited state. We also analyze the prospects for detecting gamma rays generated by 
dark matter annihilations inside Jupiter using atmospheric Cerenkov gamma-ray detectors, 
as well as the potential for observing pairs of upward going muons in neutrino telescopes 
generated by the decay of mediators to muons. Finally we offer some preliminary specula- 
tions concerning an interpretation of PAMELA'S rising positron fraction originating within 
the solar system, and discuss the obstacles to this identification as well as 'smoking gun' 
directional and temporal signatures which can confirm or rule out such an interpretation. 

The rest of this paper is organized as follows. In the next section we present the WIMP 
trapping rates and general formulae for calculating the WIMP-powered gamma ray flux, 
including angular and energy distributions. Section 3 contains estimates of the gamma 
ray signature in two variants of secluded models, with pseudoscalar('axion')- and vector- 
mediation. We conclude in Section 4 with a discussion of future prospects for improving 
sensitivity to these dark matter models via observation of gamma rays and charged particles 
in the solar system. 

2. Capture and delayed electromagnetic decays of mediators 

In this section we derive the general formulae for the expected gamma ray fluxes generated 
by metastable mediator decays. We will also consider the case of charged particle fluxes and 
comment on the expected effect of the magnetic fields of the Sun and Earth. 
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2.1 Solar capture and 7-flux 

We begin by providing the trapping efficiency of WIMPs inside the Sun, assuming elastic 
scattering of WIMPs on nuclei [HI [121 H] normalizing unknown quantities on some 
fiducial values: 

Ce 1.3 X 10- X (^52^) E/- {w^3) S('"x/™A')fA'(™x). (2) 

where the sum runs over the nuclei N present in the Sun, /jv denotes the fractional abundance 
relative to hydrogen, (^at) is the spin-dependent (spin-independent) scattering cross sec- 
tion and we have used the standard values = 0.3 GeV cm~^ and v = 270 km/s for the 
local WIMP density and velocity dispersion. The function S{x) = [A(x)^/^/(l + yl(x)^/^)]^/^, 
where A{x) = {3/2)[x/{x — l)'^]{vesc/v)'^, with Vesc — 1156 km/s an 'effective' escape velocity, 
is a kinematic suppression factor, while Fi\f{m^) determines additional suppression from a 
nuclear form-factor. This approximate formula holds for both spin-independent and spin- 
dependent scattering, but it is important to bear in mind that coherent scattering in the 
former case means that heavier nuclei, such as He, O, and Fe tend to dominate the capture 
rate despite their reduced abundance. The lack of nuclear coherence in spin-dependent scat- 
tering means that accounting for scattering off hydrogen is generally sufficient. More general 
formulae may be found in Ref. |1I] . 

For most models the trapping rate in the Sun determines the overall annihilation rate, 
as the two processes, trapping and annihilation, are usually in dynamical equilibrium. With 
an 5'-wave annihilation rate of order 1 pb as dictated by the relic abundance, a per-nucleon 
scattering cross-section of order 10^'^^cm^ or larger is generally sufficient for thermalization 
to occur over the lifetime of the Sun. Then the trapping rate Cq and the probability of a 
single annihilation event to produce gamma quanta reaching the Earth, which we denote as 
P-y, determines the overall flux of solar WlMP-generated gamma rays at the Earth's location: 



This equation is correct if the particles that produce the final photons are highly boosted 
along the radial direction of the Sun, and should be modified by another factor of 1/2 if 
those particles travel very slowly and emit photons isotropically as the Sun is effectively 
opaque and photons emitted on the back side of the Sun are reabsorbed. We are only 
interested in the former case in the current discussion. The master formula ([3]) suggests 
that observable gamma ray fluxes are indeed possible, provided that is not too small. 
In conventional neutralino-like WIMP scenarios, this probability is in fact negligible, as the 
only practical way of producing gamma rays is via high-energy neutrinos interacting with 
the outer layer of the solar materials, and in this case the probability P^ is very small. In 
contrast, secluded WIMPs offer the possibility of maximizing P^, which arises as the product 



^Earlier claims of enhanced neutralino annihilation immediately outside of the solar radius |13j were not 
confirmed by subsequent studies |14j . 
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of: (i) the probability Pout that the mediator particle decays outside the solar radius; (ii) the 
probability Br^ of producing a gamma quantum in the decay of the mediator; and (iii) the 
probability Bry of producing the mediator particle in the annihilation process: 



Here, g = g^gv where g^ and gv are the multiplicities of the photon and the mediator particle 
V being produced along the chain of reactions. If the loss of V due to re-scattering inside 
the sun can be neglected, the mediator escape probability Pout is well-approximated by 



while the branching ratios Bry and Br^ and the multiplicities are model dependent. For 
example, Br^ may be close to 1 in some secluded models, or more commonly may lie in the 
10~^ — 10~^ range for models where gammas are produced as radiation accompanying the 
decay to charged particles. 

The general relations ([2]) and use elastic cross sections normalized to 10~^^ cm^, in 
conflict with the current bounds on spin-independent elastic WIMP-nucleon scattering |15j . 
but well below the bounds on spin-dependent scattering |T7]. Therefore, the best sensitivity 
to secluded WIMPs will occur for models where nucleon scattering is predominantly spin- 
dependent. Alternatively, scattering may be predominantly inelastic xi ~^ X2, with a small 
mass gap between the two WIMP states relaxing the most stringent constraints on spin- 
independent CTp. We will consider examples in both these categories in Section 3. 

2.2 Angular and spectral distributions 

We now turn to the angular and spectral distributions of gamma quanta. We will comment 
at the end of this section on how these distributions may be generalized for charged particles, 
such as e^, for which the propagation effects of the magnetic fields of the Sun and Earth 
must be taken into account. For simplicity we shall assume a decay of the mediator into a 
pair of gammas, which are monochromatic and isotropic in the rest frame of the mediator. 
The Lorentz boost of the mediator then determines both the angular and energy resolution. 
In the absence of the boost, 7 ~ 1, one would expect the Sun to acquire a constant surface 
brightness in gamma-rays, assuming that the decays happen not far from the solar radius, 
although its side would appear brighter due to simple geometrical reasons. However, in the 
most interesting case of large boosts, 7 > 10^, the gamma quanta in the decay products are 
emitted within an angle 6* ~ I/7 from the original direction of the mediator. For 7 as large 
as 10^, the majority of gamma rays would come from an angular spot in the sky smaller 
than the solar radius, and such a tight correlation is observable with modern gamma ray 
telescopes. Also, the gamma energy spectrum would be peaked in the direction of the solar 
center, although this effect is significantly smoothed out by the finite angular resolution of 
the detectors. In this subsection, we illustrate these angular distributions and energy spectra 
for the case of = 1 TeV WIMPs annihilating into two metastable mediators, that further 
decay into 27, so that the multiplicity factors in Eq. (jlj) are gy = g^ = 2 . 



Pj = g X Pout X Br^ x Bry. 



(4) 



Pout - exp{-P0/(cr7)}, 



(5) 
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In the rest frame of the intermediate particle V, photons are emitted isotropically with 
a momentum distribution f'{p') = -;^^^5{p' — mv/2), where p = \p\. After a Lorentz trans- 
formation, the distribution function f {p) measured in the observer's frame becomes 

f{p,a) = 7(1 -/5cosa)/'[7(l - (3 cos a)p], (6) 

where (3 = v/c, and a is the angle between the momentum of the outgoing photon and that 
of the particle V measured in the observer's frame. Consequently, photons are emitted with 
an angular distribution given by l/[47r7^(l — (3 cosa)^]. The density of the photons arriving 
at the detector with an incoming angle 6 whose momenta are within a small solid angle dQ 
is calculated by integrating the particle density of V along the line of sight weighed by an 
appropriate angular distribution factor: 

d^(9) cos6'Br^ f , n(/sin^csca)/sin^ 

' dcosa 3 -. (7) 



d^l 2-KT Jd 7^ sin a(l — /3cosa)^ 

Here / ^ 1 A.U. is the Sun-detector distance and n(r) is the number density of V at distance 
r from the center of the Sun. Given our assumptions, 

^ir) = ^® exp{-r/(t;r7)} , (8) 

where v ~ c(l — 1/ (27^)) is the velocity of V . 

A few words about the range of the integral for dcosa, which we have denoted as D 
above, are in order. If the Sun were transparent, D would simply be the domain [—1, +1]. 
In the current case, however, the Sun is effectively opaque since photons traveling through 
the interior are instantaneously absorbed or degraded in energy through their interactions 
with solar material. Thus one should simply discard those photons from consideration, which 
leads to a lower bound for the integral at cos a = a/1 — {IsmOy / whenever / sin 6* < -R0. 
Consequently, 

1 - (/sin^)Vi?|, +lj , lsme<RQ] 

J (9) 
+1], /sin^>i?0. 



D 




The lower integration limits are not important for 7 ^ 1, as almost all the photons detected 
come from within a small angle ~ I/7 and the flux is completely negligible when 9 k, 9q. 
We illustrate d$/df2 for a representative value of 7 = 10^ in Fig. [2l Given a realistic angular 
resolution, it is clear that all events associated with the decays of highly boosted mediators 
will be concentrated in the angular bin covering the solar center. 

The energy distribution of $ follows from a similar integral. In general, we can write the 
full differential distribution in the form, 

d^$(^,») Z sin cos 0Br^ f , f [p, a) nil sm. 9 esc a) 

^ ^ ' dcosa . ^ ^ ^ 



p'^dpdQ 7r J^? sin^a ^-j^g^j 

lsm9cos9BT^ f dcosa \fn ,q \^ fi ■ a \ 

= / o — ■ (1 — p cosaj / [7p(l — p cosajjn^/ sm CSC a) , 

T J^i sin a 
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where D is the same domain explained above. Speciahzing to the case of two-body decays 
of V to photons, we have 



d^^{9,p) _ (72 - l)sin^cos^/ Br^ p f ^-f"^ - Isin^/p 
p^dpdn 2'K-iT{2-iPQp-pl-p^fl'^ po I ^/2-fpop -pI~p^ 



(11) 



where Pq = my/2 = m^/(27) is the photon energy measured in the rest frame of V. When 
Zsin6' > Rq, this expression is vahd as long as 



pe 




(12) 



and is understood to vanish for p outside this range. For I sin 6 < Rq, the lower limit is 



again modified to po7 ^(1 ~ /^y 1 — sin^ ^/-R© ) ^ due to absorption. 

The resulting photon spectrum varies dramatically within a tiny range of 6 that for the 
Lorentz boosts considered here is typically well below the angular sensitivity of any gamma 
ray telescope. For a detector with angular resolution poorer than the angular size of the Sun, 

the integration of (ITT]) over d^l produces a flat spectrum as long as po/'j < p < Po^j^- 
However, some detectors, such as Fermi, have an angular resolution smaller than the solar 
angular size. For demonstration purposes, we take 7 ~ 1000, and average ffTTj) assuming a 
Gaussian-profile and an angular resolution A6 of one-tenth the solar size. Since 7 is quite 
large, almost all the photons that reach the detector come from within an angle of I/7 that 
is comparable to the smallest angular size the detector can resolve. Consequently, the signal 
that would be seen by the detector represents a (relatively) bright central spot of the Sun 
with an almost exactly fiat spectrum. The photon flux drops rather quickly with angle. If 
detectable, the photon spectrum at 6* > A9 is again mostly fiat except for a very sharp peak 
near the low momentum end. We illustrate these results in Fig. [21 

Up to this point we have only discussed the spectrum of photons resulting from the 
decays of long-lived mediators. A similar analysis for charged particles, such as electrons and 
positrons, is less straightforward due to the complications of the magnetic fields of the Sun 
and Earth, as well as the solar wind, and their effects on the propagation, energy degradation, 
and absorption. A proper calculation of the angular and energy distributions is beyond the 
scope of this work, but we wish to give a qualitative discussion of these effects, paying 
particular attention to the implications regarding a possible local solar system component 
to the PAMELA signal. 

A variety of secluded models may be constructed which can lead to electrons and positrons 
in the final state. It is important to stress that the effects of the magnetic fields and solar 
atmosphere on these charged states is quite model dependent, and primarily sensitive to 
the lifetime and boost of the mediator produced in the WIMP annihilation as well as the 
production modes of the final state charged particles. If the mediators have a typical decay 
length on the order of the solar radius or slightly less, processes such as absorption and 
reflection have the potential to strongly degrade the overall signal flux, and the strong solar 
magnetic fields may drastically alter the trajectories of the charged particles. If the mediator 
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a) e b) ' P/mx c) p/m^ 



Figure 2: Angular and spectral distributions of the photon flux generated by two-photon decays of the 
mediators with jvt = 0.5Rq and 7 = 1000. (a) The normalized angular distribution ^^^^/(CQBryBr-y • 
10-2icm-2). j^^j (^j) ^Yie fully differential flux ^^^/(CoBryBr^TO-i • IQ-^icm-^) in the direction = 0, 
averaged over a Gaussian profile that mimics an angular resolution of A9 = O.16'0; (c) The same for 9 = 
0.1610. 



escapes and decays well past the solar radius, the effects of the heliosphere and the Earth's 
magnetosphere can be still be significant, especially for less energetic particles in the tens of 
GeV range and below (where the bulk of the anomalous PAMELA positrons reside). 

We must emphasize that the main feature of this signal would be a high degree of 
anisotropy in the predicted flux of electrons and positrons, correlated with the center of 
the sun, which should make this scenario testable. An important consideration is to what 
degree this tight correlation may be affected or degraded through the processes discussed 
above. It is instructive to compare this signal with a possible pulsar component to the ris- 
ing positron fraction of PAMELA [18j, where one of the main signatures may again be an 
anisotropic signal. While the latter anisotropy may be somewhat difficult to detect, this fea- 
ture in the solar component discussed here would be more pronounced due to the proximity 
of the Sun and the underlying production mechanism. This provides additional motivation 
for studies of the positional, directional, and temporal features of the existing PAMELA and 
Fermi data sets. 

In order to go further and actually attempt to fit the PAMELA and Fermi data, a better 
understanding of the predicted energy spectrum would be required, which goes beyond the 
scope of this paper. Of course, this is also quite model dependent. For example, a simple 
underlying two-body decay of the mediator to e+e" pairs will result in a very hard spectrum, 
which may not fit the observed spectrum. However, the spectrum can be softened through 
cascade decays, or the re-scattering of mediators inside the solar interior. It would also be 
important to understand to what extent the electron/positron energy spectrum is affected 
or softened by interaction with the solar atmosphere and through propagation from the sun 
to the earth. 
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2.3 Other local sources 



The Sun represents by far the largest astronomical body in the solar system, and presumably 
the most efficient WIMP capturing reservoir. However, the lifetime of the mediators can 
be such that their decays occur deep inside the solar interior and the decays outside are 
exponentially suppressed. Capture by planets may then be important and two cases to 
consider are WIMPs captured by the Earth and Jupiter. The latter option is of interest 
because of the possibility of very precise observations by means of atmospheric Cerenkov 
detectors. The flux of WIMP-generated mediators can be estimated using the basic scaling 
of ([2]) suitably rescaling the escape velocity, the mass, and the distance to the Earth. Using 
the results of [11], and assuming for a moment that V decays happen outside the solar 
interior, we find 

$j fvjY Mj (Ls-eV , ,^_9 

X X ~ 4 X 10 , (13) 



$0 V^o/ \Lj-Ej 

which represents a strong suppression. However, if the decay length of the mediator is less 
than 5% of the solar radius, the gamma ray flux from Jupiter may exceed that from the 
Sun. In this case the use of atmospheric Cerenkov detectors such as HESS, MAGIC and 
VERITAS, could set additional bounds on the capture rate by Jupiter, and consequently on 
the WIMP-nucleus cross section. 

Finally, the fluxes from the center of the Earth are difficult to present in the same compact 
form as ([2]), primarily because of the likely dependence on the annihilation rate. These fluxes 
are typically much smaller than those generated by WIMP annihilation in the solar interior 
except for some special cases with resonant energy loss [TT]. Nevertheless, the annihilation of 
secluded WIMPs inside the Earth does allow a probe of much shorter mediator lifetimes, and 
offers additional signatures in neutrino detectors. In certain models, the signature involves 
pair production of muons via the decays of the mediators, a possibility that can be efficiently 
explored with neutrino telescopes such as SuperK and IceCube. 



2.4 Observational sensitivity 

Some of the primary 7-ray observatories, such as atmospheric Cerenkov detectors are of 
limited utility in this they are unable to directly observe the Sun. While, as noted 

above, they can be used to place limits on the flux from Jupiter, a fiducial WIMP-nucleon 
cross section of 10~'^''cm^ would generate a flux that is generally too low to provide com- 
petitive sensitivity from this source. Thus, we will focus on water Cerenkov detectors and 
space-borne 7-ray observatories such as Fermi. 

The primary limit we will use here comes from the Milagro water Cerenkov detector [19] , 
which has a wide field of view and an angular resolution of 0.75°, slightly larger than the 
disc of the Sun. The background arises primarily from the scattering of high-energy cosmic 
rays, producing 7's which can arrive at the detector from within the Sun's disc. Nonetheless, 
the limits obtained by Milagro for monochromatic sources are significant and reach up to 
10 TeV. We exhibit the ensuing constraint on a contour plot of the 7-fiux in Fig. [3l which 
for characteristic scattering cross-sections, already imposes a significant constraint on the 
electromagnetic branching fractions. The primary flux limits obtained by Milagro are for 
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Figure 3: Contours (in blue) of the local 7-ray flux <i> in units of cni^s"^ in the plane of dark matter mass 
and a normalized injection cross-section assuming spin-dependent scattering. The decay distance in the 
detector frame has been fixed to Rq . The red line indicates a figure of merit for the sensitivity of Milagro 
(see the text for more details). 

monochromatic sources, while a monochromatic decay of the mediator as discussed above 
leads to a smoothed spectrum. Given that the background is a steeply-falling function of 
energy, while the expected signal is not, we adopt the Milagro bounds on monochromatic 
sources, weakening it by an order of magnitude in obtaining the bound in Fig. [31 

In the near future, the Fermi-LAT detector, which is also able to observe the Sun [20] 
with an angular resolution about 10 times better than that of Milagro, should be able to 
improve on these limits and thus provide a very significant probe of secluded dark matter 
models. The energy range of Fermi-LAT is also ideally suited to probing secluded WIMP 
annihilation with masses of a few hundred GeV, where Milagro loses its sensitivity. 

It is also possible for long-lived mediators to decay predominantly to electrons and 
positrons. It is tempting then to consider the implications of these decays in relation to 
the anomalous positrons fraction observed by PAMELA. As far as we are aware detailed 
timing and directional analyses have not yet been performed, leaving open the possibility 
of an intriguing local explanation for these anomalies within the dark matter framework. 
Given that the existing magnetic fields affect propagation, and the decay chains of the me- 
diators are model-dependent, one could in principle fit the spectral shapes of the observed 
signals. We have obtained a rough estimate of the required integrated signal flux by using 
Fermi electron data with an ~ spectrum to infer an estimate of the background and 
incorporating a new source with a harder power-law spectrum. Depending on the source as 
well as where the signal turns on (i.e. how large a component of the PAMELA signal) we 
find a required flux of 

^e+e-{E > 10 GeV) ~ 10"^ - 10^^ cm-^s"^ (14) 

We also point out that there will be an 0(1) reduction of the observed signal flux due to the 
fraction of the time the satellite is facing away from the sun. Given the assumption that at 
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least part of the excess positron flux does arise from mediator decays, this leads to a minimal 
predicted gamma flux due to the accompanying final state radiation. The typical spectrum 
of photons resulting from this process would be {a/7r)dE/E, and thus the photon flux will 
generally be no less than 0.1% of the electron and positron flux, and one infers the following 
target photon flux: 

$f > 10 GeV) > 10-^ - 10-^ cm-^s-^ (15) 

Thus, while it is possible that the hypothesis of a local origin for the flux anomalies can be 
directly tested with timing and directional information, it also appears that existing EGRET 
data [T6l [TO] may already probe the larger range of this associated gamma flux, while it is 
feasible that Fermi-LAT will be able to cover the entire range of possible fluxes, with a 
potential 10~^ cm~^s~^ level sensitivity to multi-GeV gamma rays originating from point 
sources within a year. 



3. Secluded models vs 7-rays 

Having discussed the available sensitivity to 7-rays from the Sun, in this section we will 
consider some secluded dark matter scenarios which would be subject to this indirect probe. 
Given that a large variety of model-building possibilities for WIMPs and mediators have 
been shown to exist [3], we will simply exhibit two classes of models for which the 7-ray 
flux due to WIMP trapping in the Sun may, for various parameters, be the most sensitive 
observable and source of constraints. Our approach will be to flx the parameters such that 
the decay length of the mediator is sufficient to escape the Sun, and then consider whether 
the scattering cross-section leading to capture results in a measurable (or constrained) 7-ray 
flux according to the limits discussed in the preceding section. 

3.1 Secluded WIMPs with pseudoscalar mediation 

A fermionic dark matter candidate x can be secluded by mediating its interaction with 
the SM via a pseudoscalar 'axion' fleld a with rria < m^. We will imagine a light WIMP, 
~ 10 GeV, as well as a very light mediator, rria < 10 MeV. The interactions comprise a 
series of dimension 5 operators: 

C = CsM + x(«5m7m - ^x)x + ^{d^aY - ^mla^ 

+ d^a I -^X7m75X + T^7A.75g + ^ jHl^l ] + -^aF^,F''\ (16) 

The coupling constants are above the electroweak scale but otherwise are completely 
arbitrary at this point. Since this is an effective fleld theory model, depending on the actual 
UV completion one can achieve the suppression of either f^^ or /"^ or both (see, e.g. [21]). 
The axion can be very long-lived if its mass is low enough to ensure that decays to hadrons 
and heavier leptons are kinematically forbidden. 

We will now briefly discuss the interaction rates pertinent to the 7-ray signal. We note 
that axion mediation with the WIMP sector has also been explored recently in Ref. [22]. 
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Annihilation 

Unless is parametrically larger than //, xx ~^ will be the dominant annihilation 
mode given <^ m^, and the thermally averaged rate is 



/?2 



{av) = i^-^ 2.4 X 10-26cm3s-\ (17) 



where /5 = (1 — Am^/s)^^"^ is the velocity of the WIMPs in the c.o.m. frame. Note 
that the annihilation is in the P— wave, a consequence of identical bosons with overall 
even parity in the final state. At freeze-out, taking the WIMP velocity to be 
3T|/m^ ^ 3/20, the relic density requires ~ 120 GeV x {mJlOGeVy/^. 

• Pseudoscalar Decays 

For axions in the MeV range, we can consider decays to photons and electrons, 

"^^^ " 647r3/2 ' " 27r/| ' ^ ^ 

Depending on the ratios fe/ f-y and ma/me, either the photon or electron branching 
may dominate the total width. To maximize the photon fraction, we shall assume that 
~ 10 TeV, while fe > few x 10^ TcV, in which case the decays are dominated by 
the 2-photon final states. Moreover, provided the characteristic coupling of the axions 
to quarks is small, fq > 100 TeV, the axions will be long-lived with a decay length 
sufficient to escape the Sun, 

L„ = cr„7a - 1.2 X 10^ km x (^^) ( ( V , (19) 

VlOGeV/ V rua J VlOTeVy ' ^ ^ 

where 7a ~ m^/nia ~ 2 x 10'^ for a fiducial normalization of masses and couplings. For 
the same normalization, one can explicitly check that the absorption cross section of 
energetic axions in the solar medium is too small to attenuate the flux. 

• Scattering 

The axion-like pseudoscalar mediates spin-dependent WlMP-nucleus elastic scattering, 
involving the effective axion-nucleon couplings Ca{n,p) = Wp^PJ^J^P + fn^n--)^^^n)d^a 
where l/f(n,p) = Ylq^q^'^'' / fq terms of the parameters measuring the spin-content 
of the nucleons, aI^^ = A^^^ ~ 0.8, = aI") ~ -0.5, A^^^ = A^"^ ~ -0.15. 

A straightforward calculation leads to the tree-level cross-section, conventionally re- 
expressed in terms of the 'model independent' WIMP-nucleon cross section. 



2 (1 for m„ < 4:fifjV^, 



p,n 



where fi^ is the reduced mass for the WlMP-nucleus system. Note that depending 
on the mass of the mediator, WIMP, and type of nucleus involved in the scattering. 
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there is either an enhancement or a suppression. For the trapping rate, we need 
only consider WIMP scattering with hydrogen as the scattering is spin-dependent, 
and in the specific case of the Sun, we must also account for the increase of the 
characteristic c.o.m. velocity (by a factor of 3-5) due to the WIMPs falling into the 
Sun's gravitational well. The characteristic momentum transfer in this case is |q| ~ 
2mNV ~ 8 MeV, leading us to consider very light mediators in the MeV range to 
avoid a possible velocity suppression. A light pseudoscalar with mass in the few MeV 
range will mediate an enhanced interaction due to the long range force. Effectively 
the momentum dependence cancels and we have the usual form for the spin-dependent 
cross section: 



where we have used the relic abundance to relate to m^^ and we have normalized 
the effective nucleon couplings to be consistent with the stringent constraints arising 
from rare K decays, as we will discuss below. We see that direct detection constraints 
are easily satisfied, as Eq. (1211) displays a spin-dependent WIMP-nucleon cross section 
far smaller than the (J(p^n) < 10^^^ cm^ limit [T7] . 

With these results in hand, we see that for light WIMPs and pseudoscalar mediators the 
solar 7-ray fiux can be appreciable. For such light WIMPs in the GeV range, we see from 
Eq. ([2]) that the trapping rate can be enhanced by several orders of magnitude compared to 
weak scale WIMPs, but this is compensated by a generically smaller-spin dependent cross 
section shown in Eq. fl2T|) . 

Interestingly such light mediators are not in confiict with astrophysical bounds when the 
scale of the interactions fg is below about 10^ GeV, because these states will then thermalize 
in the core of supernovae and so will not be subject to the stringent constraints from cooling. 
Constraints from BBN cannot rule out one additional thermalized degree of freedom, but in 
any case a mass in the few MeV range is sufficient to avoid these constraints entirely. Finally, 
as alluded to above, rare Kaon decays, in particular searches for vr+a, constrain fg 

to be above 100 TeV [23|, [211 [22], but pose no particular problems for the estimate of the 
spin-dependent scattering cross section in Eq. f[2Tl) . 

Using Eqs. ([21), ([3]), and ([2T|) . and assuming an C(l) branching of axions to photons 
and a decay length of order the solar radius as discussed above, we obtain the characteristic 
gamma ray fiux in the secluded model with axion mediation: 



This is a large fiux which is already close to the range probed by satellites like EGRET 
(although the mass scale is below the sensitivity range of Milagro). We have focused on light 
WIMPs and mediators, but different parameter ranges may also yield appreciable fluxes at 
the expense of some flne-tuning of couplings in the quark sector. In particular, a suppression 
of the coupling to the top quark, with larger light quark couplings, will relax the prohibitive 
Kaon decay constraints and allow an enhanced capture rate for larger WIMP and mediator 
masses. 




(21) 




(22) 
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3.2 Secluded WIMPs with vector mediation 



Unlike the axion-mediated models, secluded WIMPs lying in a hidden sector with a spon- 
taneously broken U(l)5 gauge symmetry do not require additional UV completion, and this 
sector naturally couples to the Standard Model through the kinetic mixing portal. WIMP 
scenarios in this framework are straightforwardly formulated [3], and have been a focal point 
of theoretical interest in the last year due to the positron data released by PAMELA [6]. A 
general class of WIMP models involve multi-component states Xi charged under the U(l)5' 
vector mediator, and the low energy Lagrangian after symmetry breaking involving the 
WIMP, vector V , and the Higgs h\ takes the form, 

1 1 

+ \{d.h'f-\mUh'f + '^hX + --- (23) 

If charge-conjugation symmetry is broken via /^Am, the Majorana components of the WIMP 
may be split in mass by Am ~ Amy/e', which can reduce the elastic scattering cross- 
section and ameliorate constraints on k from direct detection. The remaining particle physics 
constraints require that /t be below a few x 10^'^. 

The lifetime and decay channels for V and h' were analyzed in [23]. There are two 
regimes in which either V or h' can be very long-lived. The first refers to mh' > my and 
K < 10"^. In this case the Higgs' is short-lived, while the vector may have lifetimes in excess 
of a millisecond. This case is of no interest for us in this paper, because the trapping rate 
will scale as and will be extremely small. The second case with long-lived particles is 
rrih' < my and k > 10^^. In this case the extreme longevity of h' comes from the fact that 
its decay may only proceed at second order in the mixing angle k. This kinematic relation 
renders h' extremely long-lived even for moderately small k, while the scattering cross section 
and hence the trapping rate can remain large. The relevant interaction rates in this case are 
detailed below. 



Annihilation 

Once trapped and accumulated in the center of the Sun, WIMP annihilation may 
lead to various final states: VV , Vh', and VVV. The latter may only be possible 
when annihilation proceeds via capture into an = 1 WIMP-onium state [5], but 
otherwise this annihilation cross section is suppressed by an extra coupling constant. 
Comparison of the VV and Vh' final state branching is straightforward, once we fix 
the charge assignment for the Higgs' particle and the spin for the WIMP. If we assume 
unit charge under U(l)s for the complex Higgs' field and a fermionic WIMP, then a 
comparison of the two final states gives: 

1 < {^xx^vkA < 3, (24) 
4 (fXxx-vy v) 

where the brackets include an average over the spin orientation. The upper end of 
this ratio is achieved when the annihilation proceeds via the formation of WIMP- 
onium in which case the final state with total spin 1, decaying to Vh' is three times 
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more likely than total spin 0, that decays to 2V. The lower end corresponds to the 
case when the recombination into WIMP-onium is kinematically forbidden. Notice 
that this ratio is not changed by Coulomb (Sommerfeld) enhancement of the cross 
section, as it is identical in both channels. Thus, in this model a minimum of one 
per every 5 annihilation events results in the production of a possibly very long-lived 
h' particle that is boosted by m^/rrihr. Note that there is no strict constraint from 
ensuring the correct relic abundance in this case, as we can take this as a relation 
which fixes the U(l)' coupling a', leaving the WIMP mass as a free parameter. For 
the remainder of this section we will assume the XX ~^ mode dominates, in which 
case a' ~ 0.02 x (m^/500GeV) [3]. 

• Biggs' decays 

As noted above, the regime of interest here is when rrih' < my and the Higgs' is long- 
lived, with the dominant decay h' II occurring at order ~ x (loop factor)^ 
[2^ . Given my 3> nih' S> 2m j and a boost 7/^/ ~ m^/mh', the h' decay length is 

for mh' > 2m^ 

2 

for mh' < 2mfj_, 
(25) 

where a' has been chosen to fix the relic abundance, and we have assumed a WIMP 
mass of 500 GeV. However, the longevity of the Higgs' boson does not guarantee its 
safe passage through the interior of the sun as there is a potential loss mechanism due 
to 'inverse-Primakoff ' type conversion into V on nuclei, h' + N —>■ V + N, followed by 
prompt V decay |25]. The cross section for this process on protons can be estimated 
as follows, 

aats - 2naa'K'my' ~ 4 x lO-^^cm^ x (^^^) ' (^) ^ " ^^^^ 

This O{10~^^ cm^) scale for the cross section is sufficiently small that absorption of 
Higgs' will be negligible, but an increase by two orders of magnitude would indeed lead 
to a significant loss of h' in the solar interior. In particular, a choice of parameters 
as in the second line of fl25|) will result in attenuation of the h' fiux by more than an 
order of magnitude. Nonetheless, we see that over much of the parameter space Lh 
can naturally be large enough for /i"s to escape the Sun without scattering even for 
relatively large values of k. 

Since the dominant decays are electromagnetic, there will be significant photon pro- 
duction through various processes such as internal bremsstrahlung etc. However, it is 
of particular interest to know the branching to 7's originating from sequential decays of 
pairs of neutral pions in the product of one-loop induced h' decay, or through a two-loop 
decay directly to 77. The direct decay to 77 will have the hardest spectrum, while the 
TTo-mediated decay is likely to have larger photonic yield than internal bremsstrahlung. 
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Lh = CTh'^h' ~ lO^kmx < 



mh' / my y 

5x10-4^ UoOMeV/ 



5 X 10-3 



mh' 



100 Me V/ V 500 Me V 7 



We can estimate these branchings by constructing an effective Lagrangian as follows. 
We first integrate out all quarks that are heavier than the vector. This leads to an 
effective interaction of Euler-Heisenberg type ~ loopx (K^e^/mQ)(F^j,)^(V^^), and sim- 
ilarly a coupling between vectors and gluons. The contribution of such heavy quarks 
to the relevant decays will thus be suppressed by powers of {mv/mQY which we will 
neglect. We are left with the light quarks u, ci, s as well as perhaps c, h if they are 
lighter than the vector. The next step is to integrate out the vector from the theory 
at the loop level, which leads to the effective Lagrangian: 



Qj—hil}fil}f, 

f 



(27) 



where = —3K,'^a/2TT, and the sum runs over all quarks lighter than the vector. From 



here we can straightforwardly compute the decay h' 
the result 



77 much as in the SM, with 



I A A 'H 



647r4 



m 



V 



ml 



(28) 



where the sum is over all fermions lighter than the vector, and I is the familiar form 
factor arising from the triangle diagram. The factor in the parentheses is of order 
one. For a heavy Higgs', ruh' > 2m^ the dominant decay mode is h' — >• f^~^f^~, and 
the branching to a pair of photons is quite small, on the order of 10"'^. However if 
h' is lighter than the two muon threshold, rrih/ < 2m^ the dominant decay mode is 
h' e^e~ with a smaller total width. The branching to two photon pairs in this case 
can be sizable: 



Br 



10 



-2 



100 MeV 



for rrih' < 2,171 f^. 



(29) 



For a heavy h' it is sill possible to get a large source of gammas through intermediate 
decays h' 27r°, followed by the 7r° fragmenting to photons. The calculation parallels 
the decay of a light SM Higgs to pions, as detailed in Refs. [261 EZ] • We consider the 
theory at even lower energies, below the charm mass, where we can write the effective 
Lagrangian in terms of the trace of the QCD energy momentum tensor and match on 
to a chiral Lagrangian. It is then straightforward to calculate the partial width for 
h' 27rO: 

W..^^J(l-5^V'V(".)P. (30) 



mt 



where the function G is 





3 4^+1 
2 1 + z 



(31) 



with z = iriu/md ~ 0.56 and the sum is over the charm and bottom quarks if they 
are lighter than the vector. This function is numerically 0{\) for the parameters of 
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interest here. The 2tt^ partial width is somewhat smaller than the h' fifi mode, with 
a branching of 

Thus approximately five percent of all Higgs' decays will result in 4 photons when rrihi 
is somewhat larger than 2m^. We conclude that both for a light or heavy h' it is 
possible to obtain percent level branchings into photons. 

Scattering 

The scattering of WIMPs with nuclei in the secluded vector model was discussed in 
detail in [28], and the regime of most interest here is when a small splitting Am allows 
for Ist-order inelastic scattering on heavier elements. The interesting feature here is 
that the larger gravitational potential well of the Sun boosts the c.o.m. kinetic energy 
of the WIMP-nucleus system relative to the case for terrestrial scattering. This leads to 
a window for Am which maximizes the trapping rate in the Sun, but which is outside 
the kinematic range for terrestrial direct detection. Given ii^kin/^Af, Am/iAr -C m^: 



167iZ'^aa'K^fi% f Am 

m^ V ^kin 

which is a simple modification of scattering induced by a finite charge radius for a 
Dirac WIMP. However, this requires that Am < i^kin ~ m'^v'^ /{2fiN)- The trapping of 
inelastically scattering WIMPs, first studied in [29J was recently considered in detail 
in [301 EI] , and we can directly make use of their results in the present case. Indeed, 
the trapping in the Sun dominantly occurs through scattering on Fe nuclei, which 
are lighter than Ge used for example in CDMS. However, the larger velocity relevant 
for trapping in the Sun, means that i^kin can be larger allowing for a window in Am 
unconstrained by the direct detection limit. Indeed, the cross-section above can be 
large, 

) ~ 1 X 10-^^ cm^ X ( " \ ' ( '-^) ' (-^) , (34) 
P VSxlO-V \ mv J V500GeV/' ^ ' 

where we have exhibited the equivalent per-nucleon cross-section for scattering off Fe 
nuclei, with a' traded for the WIMP mass through the relic abundance constraint and 
m^ ^ m^r. Direct detection constraints then arise from 2nd-order elastic scattering 
which for large mediator masses in the GeV range are relatively mild, restricting n to 
be below 0.1 [28]. 

An alternative route to maximize the trapping rate while satisfying the direct detec- 
tion constraints would be to introduce several (two or more) sequentially mixed U{1) 
groups, all broken at a sub-GeV scale. This would ensure that the WIMP-nucleon cross 
section contains higher powers of and the resulting form factor suppresses coherent 
scattering on nuclei [32] ■ At the same time, due to the increased velocity inside the 
Sun, the degree of suppression in the capture rate will be much smaller. 

Given these results, the observable 7-ray flux then follows from the h' — 77 branching 
fraction Br^ ~ 10~^ discussed above and the dominant contribution from the WIMP-iron 
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inelastic scatterine cross section 10 cm to the trapping rate. This easily leads to 

a detectable flux of gammas for a wide range of parameters, 

*,e^lxlO-c.-^-x(^)^^(5^)\ ,35) 

We see that this model leads to a remarkably large gamma ray flux, that seemingly would 
have been observed by EGRET and Fermi-LAT. Moreover, for the range rrih' > 2m^, the 
accompanying muon flux will generate a flux of z/^ well in excess of the bounds set by e.g. 
Super-Kamiokande [33j. However, it is clear that generic choices of parameters may equally 
well reduce the flux (1351) by several orders of magnitude. A more comprehensive scan of the 
parameter space for this model vs the resulting flux is beyond the scope of this paper, but 
it clear that the gamma flux does indeed constitute a sensitive probe. 

4. Discussion 

Our exploratory study in this paper has suggested a number of striking indirect signatures 
associated with the annihilation of secluded dark matter trapped in the Sun (and planets), 
in the form of a high energy flux of electromagnetic particles tightly correlated with the solar 
center. In this section, we will flnish with a number of additional remarks. 

One of the issues that has become apparent is that there are relatively few competitive 
limits on gamma rays from the Sun, primarily because this is a difficult source to handle 
for many of the more sophisticated gamma-ray telescopes. However, Fermi-LAT can and 
has observed the Sun and its impressive angular resolution and long exposure times point 
to it as having the best experimental sensitivity for probing the gamma-ray signatures of 
secluded dark matter scenarios discussed here, with mediator lifetimes on the scale of the 
solar radius. Indeed, according to our analysis of model scenarios in Section 3, Fermi-LAT 
may well provide sensitivity to specific parts of the parameter space which is superior to all 
other direct or indirect probes. 

Beyond direct decays to photons, secluded models also naturally produce a significant 
branching of mediators to high-energy electron-positron pairs. As discussed earlier, its an 
intriguing possibility that massive bodies in the Solar System may provide another possible 
source of the anomalous electron and positron fiuxes in the multi-GeV range. Such scenarios 
could be probed via high-precision analyses of the spatial or temporal non-uniformities in 
these fiuxes, particularly the positron fraction observed by PAMELA. Studies of this kind 
could significantly strengthen the parameter reach in many models of secluded dark matter. 
Moreover, while the requisite fiux of charged particles could be achieved in many variants 
of secluded dark matter, they most likely will be put to the test by the upcoming Fermi 
gamma-ray data given the minimal flux that arises through flnal state radiation. 

Within this general framework, we have presented two concrete models of secluded dark 
matter consistent with the relic abundance and direct detection constraints in which the 
most accessible experimental signatures are gamma rays and charged particles from the Sun. 
The crucial ingredient is a mediator with a lifetime long enough to escape from the Sun, 
but short enough on cosmological or astrophysical scales so that there exist no additional 
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constraints arising from early cosmology beyond those already considered in the literature 
(see, e.g. [M])- It is likely that other models of this type may be constructed, and indeed 
it would be worthwhile to explore a more generic scan to understand if the gamma-ray flux 
observed in these specific models is robust. 

Beyond these points that were already touched upon in the text, we would also like to 
mention some other related issues concerning indirect detection signatures: 

• Secluded mediators with long lifetimes can, depending on the particular model and 
mass range under consideration, greatly enhance the overall neutrino flux reaching the 
Earth. For example, if the mediators are heavy enough to decay to charged mesons 
outside the solar radius, the production of muon neutrinos is inevitable. Indeed, even if 
this decay occurs inside the solar radius, the neutrino yield can be significant provided 
the decays occur outside the dense core where all muons and pions are quickly absorbed 
before they are capable of producing energetic neutrinos. 

• There is also the possibility of searching for the annihilation products of standard, 
neutralino-like WIMPs using gamma ray telescopes. If the annihilation of WIMPs in 
the solar interior creates large fluxes of neutrinos, the interaction of neutrinos in the 
upper layers of the solar atmosphere will result in the production of hadrons (ttq, K,...) 
that decay with the significant yield of 7's. The surviving 7 fraction which escapes may 
again be accessible to Fermi-LAT, and can provide alternative sensitivity to the WIMP- 
powered neutrino flux. While ground-based neutrino telescopes are ideally suited to 
exploring z/^ neutrino fluxes, the gamma ray signature will contain information about 
other flavors, and be indirectly sensitive to the neutrino energy range where ground- 
based neutrino telescopes do not have any directional sensitivity. 

• Another possibility that was not considered in this paper is WIMP annihilation to 
very light quasi-stable mediators, such as QCD axions. The conversion of axions into 
photons may occur in the magnetic field of the Sun, resulting in a 7-ray flux. This 
idea has some similarities with a recent proposal to observe the Sun's transparency to 
gamma ray sources that may contain an axionic component [25] • 

In conclusion, the existence of a generic dark sector with metastable states mediating 
the interactions between WIMP dark matter and the SM opens up new possibilities for 
indirect detection signatures in the solar system. While it is possible to search for charged 
particles produced by mediator decays, it seems the hard gamma-ray flux produced by 
decays outside the solar radius may be the most promising 'smoking gun' signature due 
to its tight correlation with the solar center, making it easily distinguishable from cosmic 
ray-induced backgrounds. As a final remark, we observe that it may be profitable to explore 
other galactic or extra-galactic implications of this scenario in which all stellar objects are 
effectively imbued with a high-energy gamma spectrum. 
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